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ABSTRACT 



The study of solar neutrinos may provide important insights into the physics 
of the central region of the Sun. Four solar neutrino experiments have con- 
firmed the solar neutrino problem but do not clearly indicate whether solar 
physics, nuclear physics, or neutrino physics have to be improved to solve it. 
Nonlinear relations among the different neutrino fluxes are imposed by two 
coupled systems of differential equations governing the internal structure and 
time evolution of the Sun. We assume that the results of the four neutrino 
experiments are correct and are concerned not with the discrepancy between 
the average rate and the predicted rate, but with a possible time dependence 
of the argon production rate as revealed in the Homestake experiment over a 
time period of 20 years. Based on the subtlety of the solar neutrino problem 
we review here qualitatively the physical laws employed for understanding 
the internal solar structure and conjecture that the interlink between specific 
nuclear reactions of the PPIII-branch of the proton-proton chain may allow 
the high-energy solar neutrino flux to vary over time. 



1 Solar Neutrino Detection: Results 



The Sun is supposed to be a simple main- sequence star that generates its 
energy through the proton-proton chain and to a much less extend through 
the CNO cycle, thereby producing a copious flux of neutrinos. The Sun is 
considered to be a simple star because very basic physical laws can be em- 
ployed to describe the evolution and the internal structure of the Sun, namely 
Newton's laws of gravity and motion, the first two laws of thermodynamics, 
Einstein's law of equivalence of mass and energy, Boyle's law and Charles' 
law of perfect gases, and Heisenberg's uncertainty principle. To probe the 
simplicity and consistency of theoretically developed standard solar models, 
over the past two and a half decades four experiments to detect the solar 
neutrino flux have been established (Bahcall and Pinsonneault 1992). 

(1) The chlorine experiment in the Homestake gold mine (USA), in opera- 
tion since 1967, observed an average 37 Ar production rate by solar neutrinos 
with energies E u > 0.8MeV of (2.2 ± 0.2)SNUQ for runs 18-109 over the time 
period 1970.8 to 1990.0 (Davis 1993). 

(2) The neutrino-electron scattering experiment in the Kamiokande mine 
(Japan), in operation since 1986 for two time periods of 1040 days (Kamiokan- 
de II: Jan87-Apr90, 450 days with E v > 9.3 MeV and 590 days with E v > 7.5 
MeV) and 220 days (Kamiokande III: Dec90-Dec91, E u > 7.5 MeV), de- 
tected solar neutrinos with event rates < $( 8 i?) b s >= [0.47 ± O.Ob(stat) ± 
0. 06 (syst)]$ ca / c ( 8 -B) of that predicted by standard solar models (Nakamura 
1993). 

For the chlorine experiment the standard solar model with the best input 
parameters predicts an event rate of (8.0 ± 3.0) SNU (Bahcall and Pinson- 
neault 1992). The experimental results referred to in (1) and (2) above re- 
vealed that the measured solar neutrino fluxes are significantly below those 
predicted by standard solar models, known as the "solar neutrino problem" . 

(3) The gallium experiment under a mountain in the North Caucasus at Bak- 
san (Russia) began to detect solar neutrinos with energies E v > 0.23 MeV 
in 1988 and reported a n Ge production rate of (851?! [stat] ± 20[syst}) SNU 
(Anosov et al. 1993). 

(4) The European gallium experiment in the Gran Sasso tunnel (Italy), in 
operation since 1991, measured a capture rate of solar neutrinos by 7l Ga 

1 1SNU (solar neutrino unit)=10 -36 captures per atom per second 



of (79 ± I0[stat] ± 6[syst]) SNU in 30 runs between May 1991 (GALLEX I: 
May91-March92) and October 1993 (GALLEX II: Aug92-Oct93)(Anselmann 
et al. 1994). 

The detected solar neutrino fluxes for the two gallium experiments almost 
agree with each other but are in conflict with standard solar models which 
predict a capture rate of (132l^)SNU (Bahcall and Pinsonneault 1992). 

The result from the neutrino-electron scattering experiment (Kamiokande) 
constrains only the high-energy 8 B neutrino flux, while the chlorine and gal- 
lium experiments (Homestake, Baksan, Gran Sasso) constrain both the 7 Be 
and 8 B solar neutrino fluxes. The results of the four neutrino experiments, 
despite of extensive experimental and theoretical efforts to reveal the origin 
of the discrepancies, do not indicate clearly whether solar physics, nuclear 
physics, or neutrino physics have to be improved to settle the account of the 
solar neutrino problem. 

The net reaction for the proton-proton chain and the CNO cycle is the 
conversion of hydrogen into helium, 

4p — ► a + 2v e + 2e + + Q,Q = 26.7MeV, 

where two neutrinos are produced in the Sun per 26.7 MeV release of nuclear 
energy. This reaction allows an estimate of the copious solar neutrino flux 
on Earth, assuming that the solar nuclear energy generation equals the Sun's 
luminosity and does not vary over time periods short in comparison to the 
nuclear time scale, 

or i 

= Q^KMMjf " 65 X 10l °^ V1 . 

where L & = 3.86 x 10 33 ergs _1 is the Suns's luminosity, AU = 1.5 x 10 13 cm 
is its average distance from Earth, and E u m 0.26 MeV is the average energy 
of the produced neutrinos. The net reaction assumes that baryon number, 
charge flavour, and energy are conserved quantities. In the following it will be 
assumed that the reported results of the four neutrino experiments are correct. 



2 Fundamental Physical Constants and Solar 
Structure: Differential Equations 

There are two distinct types of basic equations in physics: dynamical equa- 
tions exhibiting the reversible (Newtonian) time and transport equations 
reflecting the irreversible (Boltzmann-Gibbsian) time. The Sun contains a 
large number of hydrogen atoms, and evolves as a main-sequence star in nu- 
clear time scale, an extremely long time in comparison to thermal diffusion 
time (Helmholtz-Kelvin time scale) and to its fundamental pulsation mode. 
Accordingly, the time evolution of the Sun is managed by the change of 
its chemical composition governed by a system of coupled nonlinear kinetic 
equations which can be solved separately from the system of coupled partial 
differential equations of solar structure, boundary conditions, and the con- 
straint that the model luminosity at the present epoch must be equal to the 
observed solar luminosity. 

The four basic differential equations required to calculate the internal 
structure of the Sun are equations which represent respectively the distri- 
bution of mass within the Sun, the balance of gravity and pressure giving 
hydrostatic equilibrium (involving Newton's gravitational constant G), the 
outward flow of energy driven by the temperature gradient inside the Sun 
(involving the velocity of light c and Stefan's constant a), and the equation 
of nuclear energy generation within the Sun which continually replenishes 
that radiated away (involving Planck's quantum of action K). In order to 
solve these differential equations, it is necessary to specify the equation of 
state (involving Boltzmann's constant k), the opacity of the solar material 
and the nuclear energy generation rates. Boundary conditions have to be 
satisfied at the surface and at the centre of the Sun. From a more general 
point of view with regard to the calculation of the internal structure of the 
Sun it is only necessary to make an assumption about the distribution of the 
energy sources within the Sun, not necessarily knowing what produces the 
energy. Due to this fact, calculations of the internal structure of solar-type 
stars made already considerable progress before the production of energy by 
nuclear synthesis was understood in any detail (Eddington 1928, cp. Chan- 
drasekhar 1984). The differential equations for the internal structure of the 
Sun contain fundamental constants of physics that make it possible to gain 
qualitative insight into the solutions of the system of differential equations by 



simple dimensional analysis. In quantum field theory the Sommerfeld fine- 
structure constant a e i plays the role of a dimensionless coupling constant for 
the Coulomb force: 

i he 

a d l = - « 137. 

Equivalent to this quantity one can define the gravitational fine-structure 
constant a g , 

fir 

where m p is the mass of the proton. This quantity measures the smallness of 
the gravitational force between two protons, similar to a e i which measures the 
smallness of the Coulomb force between two electrons. In terms of stability 
of matter, a~l gives the maximum positive charge of the central nucleus that 
will allow a stable electron-orbit around it. It can be further shown that 
the combination of the fundamental constants in ct~i and a^ 1 are playing an 
important role in the evolution and structure of the Sun as expected from 
the differential equations (Salpeter 1966). Electromagnetic and gravitational 
interactions are long range interactions, where the unit of electromagnetic 
interaction is e and the unit of gravitational interaction is m p (m p = 1837m e ). 
The fourth fundamental constant of physics, Boltzmann's k, is an exceptional 
case since temperature can be defined in terms of energy (Gamow 1970). 
However, if the equation of state for solar material followed exclusively the 
perfect gas law, one would not have any preferred units for density and 
temperature. This leads necessarily to the consideration of radiation pressure 
(involving Planck's constant %) and electron degeneracy (involving h and the 
electron mass m e ) for the internal structure of the Sun. 



3 The Sun is Massive Because Gravity is Weak: 
Virial Theorem 

The Sun is a system of N nucleons with a mean separation of the order of 
magnitude d, 

M w Nm p , R Q « N 1/3 d, 

where M & and R & denote the mass and radius of the Sun, respectively. 
Observationally the Sun ought to be considered close to hydrostatic equi- 



librium while from the theoretical point of view the Sun is considered to 
be in complete hydrostatic equilibrium. The virial theorem implies that the 
gravitational binding energy of the Sun must be of the order of its internal 
energy, 

GM 2 „ GmlN 5 ^ 



I 



dvp = = -n a = 



Rp) d 



where P is the total pressure in volume element dV and the right-hand side 
of the equation is the total gravitational potential energy of the Sun. The 
total thermal energy content of the Sun is the kinetic energy per particle 
times its number, 

E th ~ NkT. 

Two forces balance to keep the Sun in hydrostatic equilibrium: the gravita- 
tional force directed inward and the gas and radiation pressure force directed 
outward. The total radiation energy is the product of the volume of the Sun 
and aT 4 : 

E r ~ Nd 3 aT 4 , 

where a is Stefan's constant, a = by virtue of Planck's law. According 

to Heisenberg's uncertainty principle, the Fermi momentum of free electrons 
is pf ~ h/d, where d is the average separation of electrons. Provided the 
electrons are nonrelativistic, their Fermi energy is Ep = p 2 /2m e and one 
obtains for the electron degeneracy energy, 

E H ~ N^— = N- 



2m e 2m e d 2 
Thus, the virial theorem implies that 

NkT + Nd 3 aT 4 + N-^— w N 5/3 ^^. 

2m e fl a 

For a star of solar mass with a mean molecular weight equal to 1, the radiation 
pressure at the center cannot exceed a few percent of the total pressure and 
can be neglected. Similarly, electron degeneracy does not contribute to the 
total pressure under solar conditions. Hence, 

kT ~ JV / — = N / a g — = ( — ) — , 
d 9 d \N J d 



with JV = a~ 3/2 . The Sun is massive because gravity is weak. This simple 
relationship between temperature and number of nucleons also confirms that 
the most important fact concerning a star is its mass. More detailed calcu- 
lation shows that there exists an upper limit for the mass of a stable star of 
~ 100N m p , otherwise radiation would dominate and lead to the disruption 
of it. Likewise a lower limit for a star's mass can be derived, ~ 0.1N m p , 
to account for the temperature needed to ignite nuclear fuel to form a self- 
supported shining star, i.e., to burn hydrogen (Weisskopf 1975, Carr and 
Rees 1979, Dyson 1979). 

The internal structure of a configuration M ~ lOOiVomp is dominated by 
radiation pressure, while for M > 0.1Nom p it is held together by electron 
degeneracy. 

4 Gravitationally Stabilized Solar Fusion Re- 
actor: Adjustment Factor 

The basic condition for thermonuclear reactions between charged particles is 
that their thermal energy must be large enough to penetrate the Coulomb 
repulsion between them. Nuclear reactions are collision phenomena charac- 
terized by cross sections. The cross section a of a reaction is defined as the 
probability that the reaction will occur if the incident flux consists of one 
particle and the target contains only one nucleus per unit area. The mi- 
croscopic nature of the particles requires the quantum mechanical treatment 
of the collision problem. The number of reactions is directly proportional 
to the number density of the incident flux and the number density of the 
target. In the case of the nuclear fusion plasma within the Sun, thermal 
equilibrium is commonly assumed for the ensemble of nuclei. The distri- 
bution of the relative velocities among the nuclei is Maxwell-Boltzmannian. 
The thermonuclear reaction rate is given by 

r 12 = mn 2 < (tv >i2, 

where rii and ri2 denote the number densities of particles of type 1 and 2, 
respectively, and < av >i 2 is the reaction probability in the unit volume per 
unit time. This definition of the reaction rate reveals immediately that the 
quantity 

T12 = [n 2 < crv >i 2 ] _ \ 



has the dimension of time and can be considered to be the lifetime of particle 
2 against reaction with particle 1. A suitable representation of the nuclear 
cross section contains two factors: A geometrical factor to which quantum 
mechanical interaction between two particles is always proportional, A 2 ~ 
(liv 2 )^ 1 (where A is the reduced de Broglie wave length, and fi is the reduced 
mass) and the probability for two particles of charge Z-^e and Z 2 e to penetrate 
their electrostatic repulsion: 

^ / Z,Z 2 e 2 \ 

a(v) = — -exp < — Z7r- 



hv J 

The constant S is called astrophysical cross section factor and absorbs the 
intrinsical nuclear parts of the probability for the occurrence of a nuclear 
reaction. Then, the reaction probability is defined as the product of the cross 
section o and the relative velocity v, averaged over the Maxwell-Boltzmann 
distribution of relative velocities of the reacting particles, 

f{v)dv= {^rT e ^{-£r} W ^ 

To investigate the competition between the exponential factors contained in 
the Maxwell-Boltzmann distribution function and the Gamov penetration 
factor the following order of magnitude estimation is pursued. For the num- 
ber density of the particle gas we use the mean density of the Sun with mass 
M & and radius R Q normalized to the mass of the proton, m p , 

M Q 1 

112 



The velocity of the nuclei is assumed to be the root-mean-square velocity of 
the Maxwell-Boltzmann distribution, 

(AkT\ l/2 

Vl2 = 

V m v J 

The nuclear energy generated in the Sun, which is lost by radiation, can be 
estimated in writing 

E nuc w XAmM Q c 2 , 



where X is the fraction of mass the Sun can use for nuclear energy generation, 
AmM & c 2 is the fraction of mass of the Sun really converted into radiation 
energy. Thus, the nuclear lifetime of the Sun is of the order 



E nuc XAmM Q c 2 
For the lifetime of particle 2 one has 

1 „ L & 



T 12 E; 



j rue 



Thus, 



L Q M e 2S [ 27re 2 ( m n \ V 2 



IAmM c 2 ~ i?|m pm V 2 (A;T)V2 e:rp | h G&t) j' 

and isolating the exponential term in this expression by setting it equal to 
unity and than taking the logarithm, one gets 



/ 2\ V 2 

2ira e i , « m 



2MlXAmSc 2 1 

The numerical value of the logarithmic term on the right-hand-side in this 
equation is relatively insensitive to the values inserted for the various quan- 
tities in the brackets. Using solar values for the quantities, M & m 2 x 
10 33 #,L© « 3.8Qxl0 33 ergs-\R & « 7xlO lo cm,T c0 « lO^X = 0.1, Am = 
0.007, S pp = 4 x 10~ 22 keV barn, one obtains for the logarithmic term a nu- 
merical value of about 10. Then one obtains 

kT = (^|^)-pC 2 - hkeV. 

This is the central temperature of the stationarily thermonuclear burning 
Sun. Actual central temperatures are about a factor 5 smaller or larger than 
this value due to the important fact that the majority of nuclear reactions oc- 
cures in the high-energy tail of the Maxwell-Boltzmann distribution function. 

The Sun has to adjust this temperature through the competition between the 
distribution function of relative energies of the particles and the penetration 
factor of the reacting particles. 



5 Numerical Simulation of the Five- Point Mov- 



ing Average of the Argon-Production Rate 
of the Chlorine Solar Neutrino Experiment 
: Variations Over Time 

The solar neutrino flux has been inferred from the neutrino capture rate in 
the chlorine neutrino experiment, measured over the past two decades in 86 
separate runs (1970.8-1990.0). The chlorine detector with its threshold of 0.8 
MeV is sensitive to the high energy, low flux part of the neutrino spectrum 
of the Sun, principally the neutrinos from the 7 Be and S B reactions. In the 
following we are not concerned with the discrepancy between the average rate 
and the predicted rate, but with a possible time dependence of the measured 
argon-production rate. 



Fig.l. Shows the relative importance of the different branches of the proton-proton chain: 
PPI, PPII, PPIII, producing respectively 26.23 MeV, 25.62 MeV, and 19.29 MeV. The 
differences in the energy production are due to the energy loss carried off by the neutrinos. 
The relative importance of the different branches depends on the nuclear reaction rates 
and on the temperature and density structure inside the Sun. The lifetime of the first 
particle in each reaction is indicated. 



The proton-proton chain begins by fusion of two protons (Figure 1). 
This reaction produces the great majority of solar neutrinos ($ u (pp) ~ 
10 10 z/cm~ 2 s _1 , & u (pp) ~ T^ 1 - 2 ), in which a proton decays into a neutron 
in the immediate vicinity of another proton; the two particles form a heavy 
variety of hydrogen known as deuterium, along with a positron and a neu- 
trino. The deuterium nucleus produced by the pp-reaction fuses with another 
proton to form 3 He and a gamma ray. Most often, 86% of the time, the PPI- 
branch is completed when two 3 He nuclei fuse to form an alpha particle 
and two protons, which return to the beginning of the cycle. Approximately 
14% of the time, however, 3 He instead fuses with an alpha particle, produc- 
ing 7 Be and a gamma ray; the 7 Be then captures an electron, transmutes 
into 7 Li and emits a neutrino (^ u ( 7 Be) 10 1 0vcm~ 2 s~ 1 ,& u ( 7 Be) ~ T c 8 ). 
7 Be fuses with a proton to produce two alpha particles and thus terminates 
the PPII-branch. In rare cases, about 0.02% of the time, 7 Be fuses with a 
proton to produce radioactive 8 B, which beta decays into unstable 8 Be and 
ultimately decays into two alpha particles, a positron and an energetic neu- 
trino (§ U ( 8 B) 10 7 vcm~ 2 s~ 1 , $^( 8 i?) ~ T^ 8 ). This reaction terminates the 
PPIII-branch and is extremely sensitive to the actual central temperature of 
the Sun. 



Fig. 2. Spatial distribution of neutrino sources in the gravitationally stabilized solar 
fusion reactor. The hep reaction occurs only rarely and does not influence the rates of the 
other reactions. 



Because the beta decay of 8 B follows the fusion of 7 Be and a proton, 
all nonstandard solar models predict more reduction of the 8 B neutrino flux 
than the 1 Be neutrino flux. Any reduction of the 7 Be production rate af- 
fects ultimately both the 8 B and 7 Be neutrino flux equally. This fact seems 
to make cooler Sun models incompatible with the experimental data as the 
higher Kamiokande observed rate relative to the Homestake rate cannot be 
explained because cooler Sun models reduce the expected 8 B flux more than 
the 7 Be flux. Unless there is an independent mechanism to suppress only the 
7 Be neutrino emission, all nonstandard solar models are in contradiction to 
the solar neutrino data collected in the Kamiokande and Homestake exper- 
iments. The only astrophysical explanation of this phenomenon would be a 
greater reduction of temperature in the region of the energy production (i.e. 
in the vicinity of 0.1R Q ) than just the central part where 8 B neutrinos are 
produced (cp. Figure 2). 



Fig. 3. Five-point moving average of the argon-production rate versus time (in years) 
in the chlorine solar neutrino experiment for 86 runs in the time period 1970.8 to 1990.0. 



Fig. 3 shows a five-point moving average of the argon-production rate, 
removing high frequency noise from the actual time series collected in the 
chlorine solar neutrino experiment (Haubold and Beer 1992). One notes in 
the five-point moving average that in the periods 1977 to 1980 and 1987 to 
1990 a suppression of the argon-production rate occurs. The overall shape of 
the five-point moving average suggests that there are two distinctive epochs 
spanning the time periods 1970 to 1980 and 1980 to 1990, respectively. Each 
epoch shows a shock-like rise and after that a rapid decline of the argon- 
production rate. The following simple mathematical model is able to simu- 
late just this behaviour of the solar neutrino flux based on the lifetimes of 
concerned nuclear reactions of the PPIII-branch. For the time dependence of 
each of the relevant individual reactions in the PPIII-branch, the production 
rates for argon are assumed to follow 

V = kt, 

where y denotes the argon production rate and t is the time. Then it is 
assumed that after a certain period the argon-production rate decreases to 
zero thus forming a triangle. This formulation is repeated over time. This 
type of simulation is motivated by a growth-decay mechanism discussed in 
Mathai (1993). Instead of a linear growth and decay a nonlinear growth and 
decay resulting in a bell-shaped function is also considered. That is, 

» f (t-\T + x + 2 J r\f \ 

»=-x>p( m , — }• 

where z = k = \/3 is the slope of an equilateral triangle, b = a the confidence 
level, x denotes one set of triangles delayed by x units, v is the number of 
respective triangles, and r is the half of the triangle base. The two epochs 
in Fig.3 span a time period of approximately 10 years each which will be 
divided into 40 time units of 2.5 months each. The following table contains 
the parameters chosen for the three relevant nuclear reactions in the PPIII- 
branch: For the simulation of the shape of the two epochs in Fig. 3 based 

on the parameters for the three nuclear reactions of PPIII-branch contained 
in Table 1 one has, 

/ « / (t - [q + 2iq]f \ 



nuclear reaction 
of PPIII-branch 


lifetime 


number of triangles/ 
half of the triangle 

base 


confidence level 


7 Be(e-,v) 7 Li 


10~V 


40 months / 20 


6.7 


7 Li{p, a) A He 


ktV 


8 month / 4 


1.3 


s B(e + ,u) 8 Be* 


10- 8 yr 


5 month / 2.5 


0.8 



Table 1: 



« f (t-[ r + x + 2jr]f \ 

r 5 exp i — w? — r 

s E exp \ ^— 1|- 

The result of the numerical computation for 

y(t; a, b, c; z; q, r, s; x; u, v, w) = y(t; 0.8, 1.3, 6.7; 1.7; 2.5, 4, 20; 3; 7; 4, 0) is shown 
in Fig. 4. 



Fig. 4. Numerical simulation of one 10-year epoch revealed by the time dependence 



of the five-point moving average of the argon production rate. The simulation is based 
upon the lifetimes of nuclear reactions of the PPIII-branch, producing the 7 Be and S B 
neutrinos. 
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